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Abstract 

We consider an effective Lagrangian describing a fluid living on two-dimensional planes. 
The fluid self-interacts through a Chern-Simons vector potential, whose field strength is 
proportional to the density fiuctuation. This effective Lagrangian can be related to the 
Anyon mean field, but can also be considered more generally to describe a universality 
class of superfiuids and, when charged, of superconductors. We study the relevant physical 
properties, including the spectrum, the chirality features appearing in the polarization of 
scattered EM waves, and the peculiar response under a magnetic field, i.e. a peculiar kind 
of anisotropic Meissner effect. 
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1. Introduction and summary 



We study here the physical imphcations of an effective Lagrangian a la Landau- 
Ginzburg, describing a possible peculiar superfluid. We are interested in particular in 
the case when the superfluid is charged, behaving thus as a peculiar superconductor. The 
distinctive features of this superfluid are, flrst the fact that it lives on a two dimensional 
surface, and this is in common with various theoretical studies of layered superconductors 
(systems essentially made of two dimensional layers over which the supercurrent can flow). 
Also, many high superconductors appear to behave in this way. 

The other, more peculiar, feature is chirality; that is the fact that the dynamics of 
this superfluid makes a distinction between some left and right handed behaviour. This is 
built in the deflning effective Lagrangian, which breaks parity, although the consequences 
may not be immediate to see. In fact, the effective Lagrangian contains a peculiar self- 
interaction expressed in terms of a Chern-Simons gauge field. One can thus phrase our 
investigation as a study of a universality class of some superfluids, which can be related 
to Anyon systems. Of course, one can also forget about Anyons, and consider the present 
superfluid per se as a possible physical system. 

The main point about the Chern-Simons gauge fleld is that we assume that its field 
strength (also called "curvature") is determined by the request that it is proportional to 
the density fluctuation, i. e. the actual density at a given point minus the average density. 
This assumption makes the effective Lagrangian consistent, and, in particular, compatible 
with translational invariance. 

Formally, the same effective Lagrangian has been introduced in refs. [1] [2] [3] [4], to 
provide a possible description of the Fractional Quantum Hall Effect (FQHE). Here, on the 
contrary, we do not assume a strong, uniform external magnetic field, orthogonal to the 
surface, like in the FQHE. Thus, our investigation is rather devoted to the superfluidity 
and the possible superconductivity properties. Essentially the same effective Lagrangian 
also been previously considered in ref. [5] in a generic context of Anyon superconductivity. 
Our treatment of the problem is however basically different from what discussed in ref. 
[5]. 

We devote an appendix (Appendix A) to the relation of our effective Lagrangian with 
the translationally invariant version of the Anyon Mean Field theory developed in refs. [6] 
[7] [8] [9]. 

In Section 2 we begin by discussing the effective Lagrangian and the spectrum of the 
excitations it describes, flrst by a perturbative treatment and then by studying also the 
vortex excitations. We flnd a spectrum with a flnite gap, and this gap is of course the same 
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as that found in ref. [4]. However there is an important difference in the range of values 
of a parameter, having the meaning of a selfinteraction couphng constant. Here we are 
interested in a range of values for which we show that the "roton" (or vortex) part of the 
spectrum is above the gap, contrary to the case relevant for the FQHE considered in ref. 
[4]. In fact, our range of values of the parameter is such that when the effective Lagrangian 
is compared with Anyon Mean Field (Appendix A), all the excitations correspond to the 
inter Landau levels ones, rather than the intra Landau levels as relevant for FQHE. Thus, 
in our case, the spectrum resembles the one of an ordinary superfluid, shifted up by the 
gap amount, as shown in Fig. 1. 

We then discuss in Section 3 the chirality properties. In particular we show that 
a superfluid layer has an optical activity, in that it transmits only circularly polarized 
electromagnetic waves of a definite handedness, in a range of frequencies near to the above 
mentioned gap. Qualitatively, those frequencies are in the microwave/infrared regime (for 
reasonable values of the parameters), where it would be measurably difficult to detect, and 
we show that the effect disappears at high frequencies. 

In Section 4 we discuss the magnetic properties of the chiral superfluid, discussing the 
possibility of the Meissner effect in the bulk of a sample made of a stack of many layers 
(see Fig. 2a and 2b). We find a very different behaviour, depending on the orientation of 
the magnetic field relative to the layers. If the field is orthogonal to the layers, case a, then 
the system behaves as a type II superconductor with a finite penetration length. We show 
it by considering both the penetration from an edge and also studying magnetic vortices. 
In order to get the correct result one has to include in the computation the electrostatic 
forces, which of course are present and important in a charged superfluid in general, but 
are particularly essential in our case. 

In the case 6, where the magnetic field is parallel to the layers, the behaviour is 
very different. One has to make a further distinction between to cases: the field, which 
is parallel to the layers' plane, can either try to penetrate i) in the direction parallel to 
the layers or else ii) in the direction orthogonal to the layers. In the first case i), the 
behaviour is determined, like in the layered superconductors studied in ref. [10] by the 
amount of the interlayer Josephson coupling. Namely if the Josephson coupling is small 
the penetration length in the interlayer space can be large. This is the same as for ordinary 
superconductors. 

In the second case ii) the chiral superconductor behaves differently from an ordinary 
one. We show, by an order of magnitude argument, that in this case the penetration length, 
in the direction orthogonal to the layers, is not finite and independent of the size of the 
sample. Rather, it grows with the power 1/3 of the sample's size, apart from logarithmic 
corrections. 
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Summarizing, the effective Lagrangian we have studied describes a possible physical 
system which is rich of rather interesting structures. 

2. Effective theory 

We take as a starting point the following non relativistic effective Lagrangian density 
in two space and one time dimensions: 

L = icf>*docf> - :^\D<t>\^ - 9\cf>\^ = icf>*do<t> - H (2.1) 
Am ' ' 

here (f){x, t) is a non relativistic complex field which plays the role of order parameter, 
related to the density by 

P=\4>\^ (2.2) 

We assume a fixed total number of particles A^, therefore we keep fixed N = J dPx \4>\'^- 

D is the covariant derivative D = V—iA and A is the "fiuctuation" of a CS gauge potential 
determined by: 

^ ~ 27r 

VAA^—Sp dp^ p-v^ ^\4>\'^ -v'^ v^^<p> (2.3) 
k 

As a consequence of the conservation of the number of particles we have consistently: 

J (fxV AA = (2.4) 

One can compare equation (2.3) with the standard definition of the CS field describing 

Anyons, where the CS field strength is proportional to the density. Our effective Lagrangian 

(2.1) is meant to describe the physics of the fiuctuations over the constant density, and 

therefore the dynamical effects of the CS fiuctuations as in equation (2.3). 

In Appendix A we show how to relate the Lagrangian (2.1) to the Anyon Mean Field 

theory. 

For convenience we choose to change the zero of the energy spectrum adding to the Hamil- 
tonian density the constant quantity —gv^. So that in (2.1) H becomes: 

^ = ^ I (V - + ^(101^ - = 2^ I - '^W' + (2.5) 

here we used the fact that J (Px5p = in the Hamiltonian. 



4 



2.1. Small deformations approach 



The first analysis we want to perform on Lagrangian (2.1) is based on a small deformation 
approach. Let us take the following parameterization: 

<^ = fe'' (2.6) 

Since the integral over all the plane of Sp{x,t) is zero, we can write it as the divergence of 
some quantity, and we choose: 

5p{x) = 2vV ■ u{x, t) (2.7) 
furthermore we have the gauge freedom to take u{x, t) irrotational, that is 

— * 

V A tt = dxUy — dyUx — (2.8) 
Then we have from (2.3) and (2.7): 

VAyl=^V-« (2.9) 
k 

this equation can be solved by: 

~ Anv ~ 47TV 

Ax = -^Uy Ay = (2.10) 

Notice that V ■ 1 = 0. 

With the new parameterization the Lagrangian becomes (neglecting higher orders): 



c = 2„ea„(v . s) - -v'(.veY _ _ r - -(as)^ - 4s.^(v • ur (2.11) 

Performing the variation with respect to 9 we get: 

2 2 

2^;ao(V-'u) + — V-W = =^ —Ve + 2vU=Q (2.12) 

m m 

which can be regarded as the continuity equation: 

dop + V-{pV)^0 (2.13) 
provided we identify the velocity V=^. 
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Inserting back this equation in the Lagrangian (2.11) we get: 



2m V k 



2m 



(2.14) 



so we have managed to write the Lagrangian as a kinetic part minus a potential part. 
Taking u = UQe^^^^'^^'^^ we find the spectrum of the energy: 



1 /47r'i;2\2 1 
2m 



2m \ k 



2mE^ = ^[^\ +-^(pl+plf + 8gv\pl+pl) 



(2.15) 



Here we are neglecting a contribution to the energy coming from an electrostatic interaction 
between fluctuations, which will play an essential role in the next chapter. The piece to 
be added to the Hamiltonian density is 



e 

Stt 



5p{x) [ d^x'-^^5p{x'y, 
J \x — X \ 



its contribution to the r.h.s. of equation (2.15) is +e'^v^\p\. 

The minimum of the energy, for p=0, is not zero, so we have a gap: 



E{0) = = S 

mk 



(2.16) 



Also currents can be computed: 

1 



Jqi 



2mi 
1 



(/)tD^0 - 0(D^0)t] = ^ + ^VUy^ 



^ 2mi 

which, in a compact form, can be rewritten as 



(2.17) 



Ji = -2v{ui - ScijUj) 



(2.18) 



This current exhibits the chiral property of our system. The chiral effects appear for p^O. 
To see this let us parameterize u{x, t) as follows: 



Ux{x, t) = Uqx COs{Et + ifx) Uy{x, t) = UQy COs{Et + (fiy) 

then if E = £ i.e. p — Q, the currents (2.17) can be rewritten: 

Jx = Jo cos{£t + ipo) Jy = Jo sin(£^t + (po) 
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(2.19) 



(2.20) 



where Jo and ipQ are constants depending on uq. 

Notice that this is a circularly polarized current with a definite direction of rotation. It can 
be easily seen that, since such direction depends on the sign of fc, it reverses performing a 
parity transformation. In our case k is fixed since the Lagrangian (2.1) describe a physical 
system which is not parity invariant. 



2.2. Vortex excitations 



Now we go beyond the small deformations and study the vortex solutions, i.e. vortex 
configurations, that minimize the energy: 

E = j 1^1^01^ + ^(^^-101^)^1 (2.21) 

We will follow the method of reference [11] [12] [13] based on the classical work of Bogo- 
mol'nyi [14] . Vortices are classical static solutions of the form: 

0(r, e) = /(r)e*"^ lim /(r) = v (2.22) 

r— >oo 

here n is a topological invariant corresponding to the vorticity, that is how many times the 
vortex winds round; the elementary vortex has |?T'| = 1. 

Now since in polar coordinates we have (we call here simply A what we called J[ in (2.1)): 

B = -drirAe) - deAr (2.23) 
r 

we can choose for the gauge field the form, see equation (2.3): 

t = -|i|*v(„^-/^M) 

With this choice the total energy becomes: 

^ = / + (2.25) 
from this equation we get the following requirement for the energy to be finite: 

Tl 

lim Aeir) = ^ (2.26) 

r— >oo r 



For our particular case this yields: 



27r 



drr 



—n 



(2.27) 



which is nothing but the quantization of the magnetic flux in integer factors of 2n. 

The problem now is to minimize the total energy; for this purpose it is useful the following 

identity [15]: 

\D(j)f = \{Da:± iDy)(j)f ± mV A J ± B\(j)f 



here J is the usual current 



J 



2mz 



Then we have: 



E= d^r 



|2\2 



Due to (2.3), we can rewrite the third term as: 

1 



so we are left with: 



Let us compute the contribution of the second term: 



(^^-101^)^ (2.32) 



If 



m 



1 



A J+—A = -r / de = 



2 2 

V n nv 



-n 



m r m 



Therefore the total energy is: 

,2 



E = ± n 

m 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



(2.33) 



From this expression we can learn several things. Firstly we note that the energy E is 
positive definite, see equation (2.21). We will always study the case where g>:^. 
If one wants to relate (2.1) to the Anyon mean field we notice that the value of g which 
reproduces the Anyon mean field energy is 5r = ^(1 — |) (see Appendix A) and we see that 
-(1 - l)>^ with k>2. 
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Therefore the integral on the r.h.s. of equation (2.34) is positive or zero. So we derive the 
inequahty: 

E>'^\n\>^v^ = £ (2.35) 
m mk 

So we have found that also the "big" deformations have energies over the gap £. 

Let us now study the special case g— ^ and n<0, where the bound (2.35) is saturated. 

In this case the equations simplify because the energy becomes: 

^=^'''' + 2L/^'''{'^^^"'^^^'^l'} ^^-^^^ 

which is minimal for: 

(D^ - iDy)(l) = 0. (2.37) 

Following Jackiw and Weinberg [13] we call equation (2.37) self-dual condition. 
Let us now solve the equation (2.37). In polar coordinates it is written: 

|^ + -/-^,/ = 0. (2.38) 

or r 

If we introduce the auxiliary variable 

a = -n + rA0 (2.39) 

with the properties: 



a(0) = -n > 
lim a{r) — 

r— >oo 

equation (2.38) is equivalent to the system: 

drf =-af 
r 

dra = -'^r{v^ - f) 



(2.40) 



(2.41) 



This non-linear couple of differential equations has no analytic solution; so we will study 
its asymptotic behaviours and then solve it numerically. 

As r ^ oo the equations can be linearized defining F = (v—f) — > and we get an asymptotic 
solution in terms of Bessel functions: 

F{r) = aKo( J^vr 



^ ' , (2-42) 

a(r) = a\j — rKi i.y^'"^ 
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here 



lim Ko{z) = J^ e-' K^{z) = -K'q{z) (2.43) 



and a is some constant. 
For r^O we have: 




TT 



.2„2 



n — —V r + 
k 



(2.44) 



An being some constant. 

We solved the equation numerically (with D02GAF-NAG Fortran Library Routine) and 
verified the agreement of the solution with the previous discussion. 

In the case of n>0, the equations cannot be reduced to first order and the bound cannot 
be saturated. We will study this case by a variational method which, in agreement with 
(2.35), confirms that the vortex energy for n =+|n| is in fact higher than that for n= — \n\. 
This fact should be expected since the opposite signs of n are connected by a parity 
transformation, but our Lagrangian density (2.1) is not parity invariant, so vortices that 
differ for the sign of n cannot have the same energy. 

2.2.1. Currents 

In this subsection we compute the electric current density and the total current due to 
the vortex for n < 

From (2.29) using equation (2.38) we get: 



Jr = 



Je = --fdrf 
m 



(2.45) 



The total current / passing in the plane is equal to 

ds-J= drJe{r) = / dr fdrf = / dr—f^ = 

Jo rnJc m Jq dr 



e 

m 

This current can be interpreted as the Hall current of a QHE! 
To see this fact let us consider the following Hamiltonian: 

,2 



r(oo)-r(0) =--v^ (2.46) 



H = \n\^ + ^Jd\ [\D<t>\^ + eS|0|2} (2.47) 
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where 5=— — It is straightforward to prove that this Hamiltonian is equivalent 
to our original Hamiltonian with g = ^ . In this form we can interpret the second term 
in curly brackets in (2.47) as an interaction term between the electric potential V = 
and the electric charge density e|(/)p. To such a potential will correspond the electric field 
E = —^V which in polar coordinates reads: 



Er = -—^- — — = fdrf = - — Je 

dr 2m dr mke ke^ 
r oO 



From the first of the (2.48) we read the value of the Hall conductivity cr = 



2 



(2.48) 



CT = k^. (2.49) 



2.2.2. Variational method 



We have found exact vortex-like solutions obeying the self-dual condition (2.37) for n<0. 
In order to further investigate the general case we try to find the vortices in a variational 
way. This can be done with the following ansatz for A: 



rAQ{r) = n 



1-1 + 



2 2 



(2.50) 



and minimizing the energy (2.25) with respect to the parameters p and u. Correspondently 
we get for /(r), see equation (2.22), imposing the condition /(O) = which yields the 
relation = 2a; -I- ^^^^ between the parameters, 



f\r)=v' 



— p Lor ( 
27r / 



(2.51) 



In the minimization program we have taken |n| = l, g= ^(^1 — k = 2 to compare the 
result with the minimum find above, see equation (2.35). 

2 

For n<0 we have found £^rnin = 6.44 ^ which is very close (6.44 ~27r) to the exact result, 
i.e. the gap £. 



2 

For n>0 we have found i?^i„ = 20.64 ^ which is quite above the gap. 



Notice that the variation in the number of the particles due to these vortices is: 



SN = J d\ Sp = J d\ (/2(r) - = / d^r [(^p'^r^ - 



= nk (2.52) 



So the total number of particles is preserved only if the total vorticity is zero. 
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2.3. Summary of the spectrum 



We can conclude saying that all the excitations of the system lie over the gap. 
Recalling equation (2.15), with the correction coming from the electrostatic interaction 
between excitations (see the comment after equation (2.15)), we see that for small p we 
have: 

E{p)^£+^^\p\ (2.53) 

For larger values of p it is possible to have a roton excitation [1] [2] [3] [4] which can be 
explained as a Coulomb interaction between vortices. 

In fact, the Coulomb interaction gives an additional positive contribution to the energy, 
namely the electrostatic energy due to the density fluctuation of a charged fluid. This 
electrostatic energy is obviously positive for any density fluctuation and in particular for 
the vortex-antivortex configurations. Therefore, in conclusion, the roton part of the spec- 
trum will correspond to an energy higher than the lower bound of a vortex-antivortex 
configuration, that is twice the gap £. In Fig. 1 we report the qualitative plot of the en- 
ergy spectrum versus the momentum. (In the case considered in reference [4] the vortices 
are assumed to appear in the lowest Landau level of a Hall system, whereas the gap of 
the small deformations is due to excitations to the higher Landau levels. Thus, in this 
case the vortices can have an energy less that the gap. The case studied in reference [4] 
would correspond in our formalism to ^< In our case instead we take g> see the 
discussion of section 2.2). 



Fig. 1 
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3. Chiral property 



In this section we study how the chirahty properties of our superfluid show up in 
the optical activity. Similar studies have appeared in the literature for other kinds of T- 
and P-breaking theories, see in particular ref. [16]. Let us consider an electromagnetic 
wave incoming on a thin layer of superfluid, perpendicularly to the plane of the layer. We 
begin the analysis studying how the electromagnetic potential couples with our system 
on the plane. The small deformation Lagrangian, equation (2.14), is modified adding the 

— * 

incoming electromagnetic potential A: 

1 -.2/™ 47r^;2 ^ -X" 1 2 2 2 



C = 2veV -u - —v^[Ve —u + eA] - —{Auy - 4.gv\V ■ uy (3.1) 

— * 

here u is the dual of u in the sense that Ui — eijUj. A is a two dimensional vector lying in 
the plane. All the derivative operators act here on the two dimensional coordinates on the 
plane. 

Performing the variation with respect to 6 we find: 

— = 2vV-u+—V-[Ve--^u + eA]=Q (3.2) 
68 m \ k J 

— * — * — * 

From (2.8) we have V • u = ^. Decomposing A = Al+A±, where by definition V • ^_l = 0, 
we have: 

V^ = _^^_el^ (3.3) 

V 

Then the lagrangian (3.1) can be rewritten as: 

C = -2vU ■ Ve - (V9 + cAl? - —v^ ( eA^ - ' - — (Auf- 

2m ^ ' 2m \ k J 2m^ ' 

-4gv^{V-uf (3.4) 

Now we introduce the following parameterization A = 'V(fi+'Vi(^, that is Ai = di(fi+eijdji(^. 
Thus equation (3.4)(3.3)becomes: 

• - - e'^v'^ - 1 

C = 2mif + 2veu-V^p - 2m£^u'^ + 2v8u-Vtl} - - — (V^)^ - —{Auf- 

2m 2m 

-Agv^{V ■ uf (3.5) 

Taking u = uoe'^^*+P-^^ we find: 



2 2 

uj^ -S^ - g{p)\ + 2vu ■ {iuWip + £ VV') - {Vipf (3.6) 
J 2m 



C = 2m 

here g{p) stands for the p-dependent terms in (2.15). 



By performing the functional integration over u and remembering that V<^ = ^ p and 
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y ^ _ ^ij^iPj ^ quadratic piece in A 

where: 



e// = -2^T£^ (3.7) 



^ = 2 c2 TT^ X (^1 + ^2) (3.8) 

m Lo^ — — gyp) 

The matrices Mi and M2 are given by: 

]^ ^( \iVx^+Vy^? {ip:,uj +PyS)*{ipyUJ -pxS)\ .gg. 



2„2 / ^2 



e 



M2 = ^ ( ) (3.10) 

mp2 Pi, ' 



Here we have introduced the foUowing notation 

A 



jCe// is the contribution to the total lagrangian coming from the interaction between the 
incoming wave and the two dimensional chiral fluid. Now we suppose A propagating in 
the z direction and the plane situated at 2; = 0; if we write A = A{z)e^^'^'^'^'^'^\ with A{z) 
discussed below, the total lagrangian is: 

C{A) = ^c'a^ - ^ [p^A^ + {d,Af] - ^A'fCASiz) (3.12) 
from equation (3.12) we extract the equation of motion for A: 

d'^,A+{uj'^ -p^)A-CAS{z) = (3.13) 

For A{z) we have: 



= I (3.14) 



Imposing continuity at ^ = for A{z) and its derivative, we get: 

a+ = -a-, P- = q;+ — q;_ = -a- (d.lo) 

2ik-C 2ik-C 

Prom equation (3.15) we can read the matrices of reflection and transmission: 

1 

R^—^ik^ T= ^ (3.16) 

1 — — 1 — — 

^ 2ik ^ 2ik 
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Notice that R = T-1 and |Tp + |-Rp = l as it should. 

In particular we are interested in the coefficient of transmission in the limit when a;— >f . 
In this case only Mi in equation (3.12) is important. If we compute the inverse of the 
matrix 1 — 2^ £ we get, taking for simplicity p= (0,p): 

^= 1.2M,l,^.2,.s {'7 l-i(7+2.))+«(V.) (3.17) 

here we have set q = -^^^^ and s = 
The limit u^S corresponds to qf— >oo that is: 

T=It^(] (3.18) 



2i + s \i 1 
Now let us consider the following parameterization: 

Aj; = i?e(cos^e~*'^*) = cos^cosa;* 

Ay = Re{sinee''^e-''^*) = sin^cos((^ - ut) 

we can see that the circular polarization is for 0=j and (fi=^, so that: 



(3.19) 



V2 . V2 , , 

Ax = -— cos ujt Ay = -—siacvt (3.20) 



in our complex formalism this corresponds to: 



"''nL)e-'-' = 4(^-]e-'-' (3.21) 



Thus we see that T projects a state in a circularly polarized one, that is we can write: 

1 i 

2 i + s 



T=--^\c){c\ (3.22) 



where 

\c)=(]] (3.23) 



That is, at the resonance uj^S our chiral planar system behaves like a perfect polarizer. 
Finally note that if a;— > 00 we get: 



e^v"^ (I 



SO we have lost any polarizing effect. 
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If conversely we take a;— >0 and p^O we get C = 0, so there is no more coupling between 
the electromagnetic wave and the planar system. 

4. Coupling with an external magnetic field 



4.1. Magnetic field orthogonal to the layers 



4.1.1. Meissner effect 

In this section we test our chiral charged superfluid with an external magnetic field, and 
study its Meissner effect. 

Up to now we have been considering a two dimensional system; now we are going to 
study an effect which is essentially three-dimensional. Therefore we suppose to have a 
multilayered bulk of many two-dimensional thin films separated by a spacing d (for a 
review on the properties of the layered superconductors see, for example [17]). We further 
suppose that at the edge of the bulk there is a uniform, constant magnetic field orthogonal 
to the layers' plane as in Fig. 2a. 



Fig. 2 
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In this three dimensional system we must take into account also an electrostatic contribu- 
tion, (see the comment after equation (2.15)) and which will play an essential role here. 

— * 

This contribution is essentially due to an electric field E which comes from the fluctuations 
of the charged matter and obeys the Maxwell equation: 

V-E = e5p^^^ (4.1) 

here Sp^^^ = ^, is the three dimensional density. 

In other words if there is a fluctuation Sp of matter the system will not be in electrostatic 
equilibrium anymore, for there will be some zones where there is lack of charged matter and 
others where there is abundance: E is the electric fleld resulting from this non-equilibrium 
situation. 

So let us consider the following three-dimensional Hamiltonian density (in this section we 
have redeflned A^A'-''^ to distinguish from A'^'^): 

here 

S = VAA'='" i = -VAg'" D = V -il^^ - iel^"^ (4.3) 



Notice that the following equations hold: 



V • ^ = |5p f V • A^^ = 



(4.4) 



here ^ is a three-dimensional vector whereas A^^ is two-dimensional. If we suppose that 
the matter distribution is constant in the z direction, which means that we have exactly 
the same matter distribution in every layer, then we have no electric fleld orthogonal to 
the x-y layers' plane. With this assumption equations (4.4) tell that E and A'~^^ are dual 
two-dimensional vectors, i.e. : 



ke 

r 

So we can write the hamiltonian as follows 



E. = ^e,AOS (4.5) 



^ = [ dxdyl-^\Dcpf + l 



Lz .1 I 2md 2 



+ ^lj,-vY't (4.6) 



here we have performed the integration in z. 

At this point we make one further assumption: we suppose that the external magnetic 
fleld is constant in y, which is the direction parallel to the edge of the bulk, so we can take 
A^'^{x) in the y direction: 

A""^ = {0,A''"'{x),0) (4.7) 
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furthermore we still have the freedom to take also A'~^^{x) in the y direction: 



A^^ = {Q,A^^{x),Q) (4.8) 

In the gauge: 



y.^-em^^.^CS^Q (4.9) 

j^em g^j^j j^cs depend only on x. With these assumptions all quantities in (4.6) depend 
only on a;, so we can perform the integration in y and get: 

{A'^'f + ^-ip-v^f^ (4.10) 



+ 



Ly being the length of the edge of the bulk. 

Now, before going into computational details, wc want to spend some time analyzing 
qualitatively the reason why we expect a Meissner effect. To this end we will make some 
simplificatory assumptions. The first is to suppose that some external device is keeping 
constant the magnetic fliix $o — Lyl^B, where Ly is the length of the edge of the sample 
and Ix is the penetration length of the magnetic field, so that 

Bla^ = ^= ipo (4.11) 

Ly 

is given to the system from outside. 

The second assumption we make is to suppose B constant for < a; < Z^; and zero for x>lx 
whereas actually it is exponentially decreasing. 

The third assumption we make is to take a constant value for the density p = v^. 
So we have for 0<a;<Za;: 

B=^ A<^^ix) = f (^0 (4.12) 

''X ''X 

For x>lx, A^"^ is constant and it is possible to cancel it in the covariant y-derivative with 
the phase of (p. Thus, the integration runs from to Z^;. 

To test the meaningfulness of what we are doing let us see what happens for the well known 
case of the standard superconductor, i.e. with A'~^^ = 0; the hamiltonian (4.10) becomes: 
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Minimizing H with respect to Ix we find Ix = \J~^^ which is indeed of the same order as 
the standard value of the penetration length for type II superconductors: 



Notice that in what we have done a fundamental role is played by the term quadratic in 
^•""j the "mass term" for the electromagnetic field. 

In our case, with A^^ ^ 0, this effect could be ruined by the possible cancellation A'^^ = 
_g^em 1^^^^ notice, we have also the electrostatic interaction term quadratic in A'-^^ which 
now plays the dominant role. We have: 

Minimizing with respect to Ix we get Ix = ^^r^ we see that we get a finite penetration 
length, so we expect to have Meissner eff'ect. 

Now we turn back to hamiltonian (4.10) and analyze it more quantitatively in the 
framework of the small deformation approach introduced in the previous chapter. 
So we have, recalling the basic definitions, see equations (2.6), (2.7) and (2.10): 

(j) = f dx^{x) = dlu{x) (4.16) 

5p = 2vdxu{x) (4.17) 

Notice that in this one-dimensional case (f) is real, so there is no phase. 
The hamiltonian becomes: 



H 1 ( 1 

ax 



LyLz d Jq 1 2m 



+ 



40^A'- f + H^^n^ + 1^ f 1 - i ) O^nf \ (4.19) 



2 d m \ k 

We now make the following ansatz: 

u{x) = uoe-^'' A'^'^ix) = Aoe-^'' (4.20) 
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here uq and are the values at the edge of the bulk. The resulting equations for u{x) 
and A'^'^{x) can be written in the form: 



mk 



mk 
2m 2^ 



Putting to zero the determinant of this matrix we get: 



d 

Am' 



2nv^d / 1 
m \ k 



+ 



4m? 



4 4 

e V 
md 



Snv^d 
2mk^ 

= 



u 

J^em 



m^ 



= 



(4.21) 



k 



(4.22) 



If we take the following typical values for the parameters (as an order of magnitude 
we take m to be the electron mass): 



d=lk 



137 



= 4 • lO-^A-^ 



m 



= 250A-^ k = 2 



we get the solutions: 



Ai = 1.21 • 10 



-3 



A2,3 = 0.56 ± iO.54 = a±ip 



(4.23) 



(4.24) 



We can approximately compute analytically the smallest eigenvalue Ai by rewriting equa- 
tion (4.22) in an approximate form taking only the leading terms: 



A^ 



_d_ 

4m 



-A^ 



2nv^d 



m 



\2 I 2 2 

A + e V 



4 4 

e V 



md 



(4.25) 



for small A^ we get the solution: 



X' = 



e V 



md e^v^ md 



(4.26) 



which, compared with equation (4.14), is exactly the expression for the inverse of the 
square of the standard penetration length, and, corresponds numerically to Ai in equation 
(4.24). 

Let us note that by solving (4.22) in the limit of very small we would find: 



A 



4:TTd 



(l + 0(e^)) 



(4.27) 



which, apart for inessential numerical factors, is equal to the inverse of the Ix found mini- 
mizing (4.15). Indeed the qualitative analysis leading to (4.15) was meant in the limit of 
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a very small electrostatic interaction. Actually this is not the case. For instance, with the 
values given in (4.23) we get ^^<^e\^, and (4.26) follows. 
The general solution for u{x) is the linear combination: 

u{x) = «ie-^i^ + uae"^'"^ + use-^^^ (4.28) 

Imposing the reality condition u*{x) = u{x) we can write: 

u{x) = ttie"'^^'" - e~°''^{ui cos^x -wsin^x) (4.29) 

giving for 6p: 

Sp = 2vdu = 2v -AiMie-^i^ + e"""^ {{aui + (3w) cos l3x + {l3ui - aw) sin l3x) (4.30) 
Requiring p to vanish at the edge, i.e. Sp{0) = —v^ we can determine the value of w: 

The following relation holds between A^'^{x) and u{x): 

A'^{x) = ^^—^-^u{x) (4.32) 
^ ^ mkd A2 - ^ ^ ^ ^ ^ 

ma 

We then compute B = dxA^^. 

Imposing -B(O) = -Bo we get a value for ui, the only parameter still undetermined. If we 
substitute the values of the parameters we see that the exponential behaviour of B is 
controlled by the first eigenvalue, i.e. the penetration length is ^ = 826. 45A. We can 
compare this value to the numerical value of the standard penetration length (4.14) which 
is 825.57A. 

It is seen numerically that the behaviour of 5p is controlled by the other eigenvalues, i.e. 
the coherence length is ^ = 1.78 A. 

Notice that the penetration length is about 400 times the coherence length, so our system 
behaves like a type II superconductor. 



4.1.2. Vortices 

In this section we study, with a variational method, the vortices in presence of an external 
electromagnetic field orthogonal to the layers' plane. 

These vortices have an origin different from the vortices studied in section 1.2. Those were 
originated by fluctuations of the CS magnetic field, i.e. of the matter density from the 



21 



mean value. The present ones are instead the standard well known vortex configurations 
of superconductors of type II between the two critical temperatures. They consist of 
small regions of the specimen of normal behaviour, all surrounded by a superconducting 
region, where the external magnetic field penetrates completely and uniformly. There are 
superconducting currents fiowing around the vortices. Furthermore there is a penetration 
of the magnetic field from the vortex region to the surrounding superconducting region 
(Meissner effect). 

As we saw in the previous section the hamiltonian is: 

where we have included the electrostatic effect, that is the third term on the r.h.s., as in 
equation (4.6). 

We make the following ansatz: 
= /(r-)e^"' 

r^c-s = ^ r dr'r'[f{r')-v^] = ^^V^e"' 
k Jq k 

from the last one we easily get: 

/2(r) ^v^-v^{l- a;r2)e-'^'"' (4.35) 

Notice that our ansatz is such that we have quantization of the flux of the external mag- 
netic field, differently from the vortices studied previously when it was the flux of the CS 
magnetic field to be quantized. Notice also that since the vorticity is no longer connected 
with the CS field now, differently from the vortices studied before, it is possible to have 
an isolated vortex or antivortex since drr[v'^ — /^(r)] = as can be checked from the 
last of equations (4.34). 

If one substitutes in (4.33) and minimizes, numerically, with respect to A and u, with the 
usual values of the parameters (4.23) one gets: 

CO = O.lOlA-2 n = 7.336 • 10"^A-2 (4.36) 

for the case of the vortex (n=l) and 

oj = O.lirA-^ II = 7.337 • IQ-'^A-^ (4.37) 



(4.34) 
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for the case of the antivortex (n=— 1). 

From these results we get the values of the dimensions of the vortex and of the fluctuations 
of the magnetic field: 



1 



= 3.147A 



= 1167.54A 



(4.38) 




for the vortex, and 



1 



= 2.924A 



= 1167.46A 



(4.39) 




for the antivortex. Notice that for the vortex and for the antivortex we have difi'erent 
results, although very close to each other, as it should have been expected since our system 
is chiral. 

We see that again the behaviour is that typical for a type II superconductor. 



We can also study the properties of these vortices due to the fact that the superconducting 
material is composed of a stack of many not-strongly-coupled layers. We follow reference 
[18] where a three dimensional vortex is built up superposing a stack of two dimensional 
vortices in the case of a standard high superconductor. In that paper it is explicitly 

— * 

computed, as a first step, the magnetic field, 6, produced by a single layer and then, using 
this result, the whole stack contribution is computed. Here we will only show that in our 
case we can reach, under reasonable assumptions, the same first step, and then just state 
the final results. 

The problem has a cylindrical symmetry therefore we will use cylindrical polar coordinates 
(r, ^, z). The supercurrent flowing in our single layer is: 



4.1.3. Layered structure of the Vortices 



J=(0,Je,0) 



Je{r,z)^Ke{r)b{z) 



(4.40) 



From Ampere's law we get: 




(4.41) 



from the last two we get 651 = 0. 

For equation V A 6 = is solved by: 




— Ag — 







(4.42) 
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notice that V • ^ = 0. Ao{q) is to be determined. 

Now from the first of (4.41), solved for z = 0, and (4.42) we can get: 



Ke{r) = br{r, 0+) - br{r, 0") = 2 / dqqAo{q)Ji{rq) (4.43) 

^0 



For our vortex we know: 



e /n 



Ke^-r-- eAr - A^^ ) (r ) (4.44) 



therefore, using (4.42): 



^0 



dqAo{q)Ji{rq) 



2q+^fir) 
m 



e 

m 



Now if the superconductor is of type II i. e. if the dimensions of the vortex are much smaller 
than the penetration length of b, then it is sensible to take /^(r) equal to its mean value 
v'^ and ^^'^ = 0. Within this approximation (4.45) becomes: 

r dqAoiq)Mrq)(2q+^) = ^(^) (4.46) 
Jo \ m J m \r J 



Using the orthogonality property of the Bessel functions drrJi{rq)Ji{rq') = ^5{q — q') 
we get: 

AM = 2g+ 4.47 

A is related to the nominal penetration length = J^(3) by the relation: 

A = ^ (4.48) 

and can be considered as the two dimensional penetration length. 

One can then compute (4.42) and verify that the magnetic field decays with the penetration 
length A. 

These results are exactly the same obtained in [18]. We will not reproduce here all the 

computations that can be found in [18], but just state the main results. 

Having studied what happens with a single layer we possess the building block to all the 

multilayered system superposing the entire stack of two dimensional vortices. 

Then one can study the binding energy between vortices in different layers and find that 

thermal excitation breaks up the stack above a transition temperature corresponding to 

the Kosterlitz-Thouless temperature for a bidimensional system, see [19] [20]. 
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4.2. Magnetic field parallel to the layers 



4.2.1. Meissner effect 



In this section we consider a different physical problem with B in the plane of the layers. 
This problem, for the case of standard high superconductors, has been studied in many 
paper by J.R. Clem and collaborators [10] [21] [22] [23]. In our analysis we choose B point- 
ing in the x-direction (see Fig. 2b), and study the penetration length in z, Ig, supposing 
uniformity along y. Let us start from the following Hamiltonian: 



L, J 



dxdz 



1 



2md 



\lB'+9{p-vr}- 



— / dxdx'dzdz' 
2 



{x — x'Y + — 5^(2^5 z 



5p{x' , z') 
d 



log 



A2 



d 



(4.49) 



here we have introduced the electrostatic interaction between the fluctuations of the 
charged matter which, due to uniformity in the y-direction, is in fact the electrostatic 
interaction of a two-dimensional distribution of charge, which is known to be logarithmic; 
A is a convenient dimensional constant which, since J 5p = 0, can take an arbitrary value. 
In the gauge where A^"^ and A'^^ are both in the y-direction — (0, A^'^{x, z),0), A'-'^ = 
(0, A'-^'^{x, z), 0) we can rewrite the Hamiltonian as: 



—— = / dxdz 



+g{p - v'f 



2md 

2d2 



+ 



md^jo 



1 



cs 



dxdx' dzdz' < 5p{x' , z') log 



2md 

{x - x'Y + (2; - z'Y 

A2 



bp{x,z) \ (4.50) 



here jo is a constant depending on the intrinsic features of the material and measures the 
Josephson coupling between neighbouring layers. 

Similarly to what we have done for the field orthogonal to the layers we consider the 
possible cancellation AP^ix^z) = —A'^'^{z)[9{x)9{La; — x)], being the length of the 
x-edge of our sample, and 9 is the step function. Thus: 



Sp{x, z) = ^d^A^'ix, z) = -^A-'^iz) [5{x) - 5{L, - x)] 



(4.51) 



The first and the last term in the first integral of (4.50) are edge effects giving a contribution 
which is very small compared with the rest since, going through the computation, one can 
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see that they are suppressed respectively by factors (L^) ^ and (Lx) So in the following 
we will skip them. Thus: 



1 + 



(z - z') 



l\2 



A^'^iz') )■ (4.52) 



Like we have done in the previous section (compare with equation (4.15)), we fix the total 
flux ^Q = BLylg and define (fQ = such that for Q<z<lz'- 



B 



(4.53) 



and the integration is from to Z^. 
Therefore we get: 



fc2e4 1 



dzdz' < zz' log 



1 + 



{z-z')^ 



+ 



2L IdTT^d^ 



(4.54) 



We assumed Lx/lz large to expand the logarithm. 
Minimizing H with respect to we get: 



37r- 



2 \ 1/3 



1/3 



.og( ^ ) + 



n-1/3 



(4.55) 



oo 



Therefore we have a "quasi Meissner effect" in the sense that even if — > oo for 
we have — 0. The result of equation (4.55) appears to be quite peculiar of the chiral 
superconductor studied in this paper. 

In fact, the non-chiral superconductors would show ordinary Meissner effect for the ge- 
ometry considered in this section (at least for not too small Josephson coupling between 
neighbouring layers): that is, magnetic fields parallel to the layers would penetrate in the 
direction orthogonal to the layers (say the ^-direction) for a finite distance Iz, of the order 
of the standard penetration length. 

The chiral superconductor studied in this paper would instead penetrate much more, for a 
distance Iz as given in equation (4.55). This result holds for arbitrary Josephson coupling 
between neighbouring layers. 

(Let us also mention that, in the case of very small Josephson coupling it is possible to 
imagine another configuration possibly giving a large penetration in z. This configuration 
is not peculiar of the chiral superconductor studied here. We discuss this configuration in 
Appendix B). 
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Finally we review a general argument for multilayer superconductors indicating that the 
magnetic field parallel to the layers can penetrate easily in the interlayer spacing, therefore 
in the direction parallel to the layers (say in the y-direction, for a field along x) , even in the 
case in which Iz is finite. How much the corresponding ly is large, depends on how small is 
the interlayer Josephson coupling. This effect would be basically the same for chiral and 
non-chiral superconductors. In order to study this penetration in the direction parallel to 
the layers we follow the study of ref. [10], considering a possible vortex configuration inside 
the material. From this configuration we will infer the penetration properties. A crucial 
fact for this study is taking into account the Josephson coupling between neighbouring 
layers [10] [24]. This study has already been worked out in [10] for the case of a standard 
high Tc superconductor, here we merely rephrase that paper for our case. The main fact 
is that one has to add to the supercurrent density: 

2 

J=— (V7-eyl) (4.56) 

m ^ ^ 

(here 7 is the phase of (p, equal to in9 for the vortex configuration) , the Josephson current 
fiowing between to neighbouring layers, say the n-th and the (n + l)-th, proportional to 
the sine of the gauge invariant difference of phase of between the layers: 

rn+l 

jz = jo sin(A7n) A^n = ln+1 -In + e dl-A (4.57) 

J n 

Jo is a constant depending on the material. 
Solving the Ampere equation one finds [10]: 

with Aj = here d is the stack periodicity of the layers. 

This result tells essentially that B penetrates differently along y, i.e. parallel to the layers, 
and along z, i.e. orthogonally, in other words the vortex has an "elliptic symmetry" . Notice 
that in reference [10] the decay length along the y-axis is called instead of Aj; the reason 
is that this decay length is essentially due to the screening currents which are pointing in 
the orthogonal direction i.e. the 2;-direction. Similar argument goes for the decay length 
along the 2;-axis which in [10] is called A^. 

In conclusion, in both cases of ordinary and chiral superconductors the magnetic field 
would easily penetrate along the interlayer spacing. Instead, the penetration across the 
layers of a magnetic field parallel to them will be finite for an ordinary superconductor 
(for not too small Josephson coupling), whereas it will grow with a fractional power of the 
sample's dimension for the chiral superconductor. 
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Appendix A. 



In the standard Chcrn-Simons (CS) description of Anyons one introduces a CS gauge 
field A, whose field strength is proportional to the density: 

VA^'ocp (A.l) 

In this formulation to have translational invariance one has to take care of the boundary 
conditions in a proper way and to take a torus [6] [7] [8] [9]. In this way one shows 
that the ground state of the full quantum solution of the mean field theory corresponds 
to a constant density. Due to the non-trivial topological properties of the torus one has 
to take into account the topological components of the gauge potential (also called "fiat 
connections" since V A a = 0) defined by: 



o-x = ^ dxAx ay = <j> 



dyAy (A.2) 



here the two integrals are performed along the two non-trivial loops of the torus. 
In this way the Hamiltonian can be written as follows [6] [7] [8] [9] : 

where the covariant derivative is: 

Di ^ di- ^'^'^^ij^j ~^J^ ~ ('^•^) 

here L is the length of a side of the torus (we are supposing, for simplicity that our torus 
is a square with identified edges), is the mean density, A is the fiuctuation part of A 

— * 

— * 

such that V A^oc^p. c is a positive constant whose value can be arbitrary in what follows. 
In the Anyon problem the value of c is large, c— > oo, and thus the degrees of freedom 
remain in the ground state. 

In the case of the infinite plane we have to take the limit L — > oo at constant density (we 
will call it thermodynamical limit), so ^ ^0 if is bounded (we will see below that indeed 
tti is bounded) but in the Hamiltonian still survives a term in a^. 
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Note that now the wavefunction is now a function of a beside x; i.e. '^ = '^{x,a) so we 
define the density p{x) to be: 

p{x) = J d^a\iP{x,a)\^ (A.5) 

It is convenient to introduce a complex notation: 

A = Am -iA A = Am + iA (A.6) 



where 



are the "mean field" parts and 

A='-{A^-iA2) A = -'-{A,+tA2) (A.8) 
are the "fluctuating" parts; for the flat connections we define: 

a=7^(iai + a2) o = 77-(— ioi + 02) (A.9) 

ZTT ZTT 

For the CS magnetic field we have: 

27r 

B = 2i(dA-dA)=BM + B = -—p (A.IO) 

k 

where: 

Bm = -fv' B = {dA- dA) = '^iv^-p) = -'^Sp. (A.ll) 
Let us now suppose that ^4 = and study the mean field solution. The Hamiltonian is: 

The state of minimal energy corresponds to: 

V' = V'M = e-^^^"-^"«^(^,a) (A.13) 



here g is an arbitrary antiholomorphic function. 
From (A.13) we can see that a is bounded. 

Let us look for the constant density solutions, which are known to correspond to the ground 
state of the full quantum mechanical mean field problem. Choosing g{z.,a) = e^'"^"' v we 
get Pm = v'^, in fact 

PM = j (fa\^M\'' = j (fae--^\''-l'\\^ =v^ (A.14) 
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provided we normalize the measure (Pa such that J dPae = 1. 
For the Hamiltonian (A. 12) we get: 



TT 



Hm = —rV'N 

mk 



(A.15) 



this is the energy of TV particles in the lowest Landau level. We know that in the mean 
field solution actually the N particles fill exactly k levels corresponding to the energy: 



Em = ^\BM\Nk 
2m 

So, in order to reproduce the correct mean field energy, we must add: 

B' = _L|B„|iV(.-l) = |(l-i)/A/d^.W 

So our correct starting Hamiltonian is: 



(A.16) 



(A.17) 



(A.18) 



Notice that more in general a constant density is also obtained taking: 



M 



e 2fc 



■ zz — ^aa+nvza+i{pz+pz) 



then the Hamiltonian (A.18) becomes: 



Hm = d z< — \ip\ V 
m 



(A. 19) 



SO we have found, beside the standard mean field energy (A.16) a kinetic energy equal to 
that of one particle times N. So our system moves like a condensate where all particles 
have the same momentum. In other words it represents a collective motion. If we compute 
the currents: 



J 



-! 



1 



2mi 



we get: 



J 



m 



-p 



(A.21) 
(A.22) 



So we have found for the current exactly the charge density times the velocity. 



Now let us introduce the fluctuations taking A^O: 



H = 



cfz 



d a i — 
m 



c 

+ — 



d 7rk_\ , 2 



(A.23) 
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Then if we take: 

for the density we get: 

J £a\iP\^^ J (fae-''^\~^-^^\^\4>{z,z)\^ = \(t){z,z)\'' (A.25) 
and for the Hamiltonian: 

H = I rfV|^|(V-a)0|V^|0r} (A.26) 

where in order to reproduce the correct "mean field" energy we have to take g = -^{l — 

So we have found an Hamiltonian which in the "mean field" case gives the correct 
answer and in the more general "fluctuating" case recovers the covariant derivatives taking 
into account the density fluctuations. 

We take equation (A.26) as the basic effective Hamiltonian, describing our quantum fluid, 
corresponding to the effective Lagrangian density (2.1). 
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Appendix B. 



We consider the geometry of section 4.2, Fig. 26), and a configuration relevant for the 
penetration of the magnetic field in the 2;-direction. 

Assume that the system is homogeneous in the a;-direction. We suppose that (j) is constant 
up to a phase i.e. (f) = ve^'^, therefore we have no electrostatic interaction and put ^'-^'^ = 
(therefore this configuration would be there also in a non chiral superconductor). 
In this configuration one considers the possibility that the phase is responsible for the 
cancellation of A'^'^: 

dy^ = eA^"^ (B.l) 

We again suppose to have a fixed external flux so that equation (4.53) still holds. 
We further impose periodic boundary conditions (ensuring the current conservation): 

7(0) = 7(Ly) = O(mod 2tt) (B.2) 

From (B.l) and (B.2) we get: 

eLyA^"^ = 27TP (B.3) 

here p e Z has in general a different value in different layers. Since d is the spacing between 
the layers the z coordinate of the s-th layer is Zs = sd. Now recalling equation (4.53) we 
get: 

eLyA'^'^izs) = -eLy^sd = 2ttPs Ps = -^^(foLyS e Z (B.4) 

iz 27r tz 

but since the total flux is quantized (fioLy = we have: 

— = --n (B.5) 

S Iz 

since ^ is the number of layers "penetrated" by B this means that every layer bears an 

identical fraction of the total flux. 

The Hamiltonian of this system can be written as: 

H^J dxdydz!^^\dy^-eA^-\\^Jd,^\' + l{dzA^-^f^ (B.6) 



32 



the second term is the Josephson coupling between neighbouring layers [24]. Aj is a 
constant of the dimension of a length whose physical significance will be clear in the 
following, it is related to the Josephson maximal current density jo by the equation: 



Jo 



(B.7) 



A is the standard penetration length (4.14). Since from (B.l) we have j = eyA'^'^ , we can 
rewrite the hamiltonian as: 



1 /L„ 
1 + t: 



3 V A 



2n 



(B.8) 



which is minimal for Iz macroscopic i.e. for Iz ^ Lz ^ Ly L^- To see if this configuration 
is favoured one has to compare it with the energy of the standard case in which: 



H 



standard 



^0 



2 LxLy 



(B.9) 



That is, for Iz^Ly, the configuration giving (B.8) is favourite if: 



1 1 
TIa ^ 2/1 



1 fLy 

1 + 



3 V A 



(B.IO) 



For all common high materials Aj ~ lO^A (corresponding to j'o — lO'^A/cm^, see for 
example [25]) therefore since ^ ~ 120, for Iz ~ lO^A, this configuration could be 
marginally competitive. When equation (B.IO) is not satisfied then this configuration 
would not be favourite. 
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Figure Captions 



1. Qualitative spectrum. 

2. Directions of the magnetic field with respect to the layers. 
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